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Abstract 

We study thermodynamics of A/" = 6 superconformal Chern-Simons theory by 
computing quantum corrections to the free energy. We find that in weakly 
coupled ABJM theory on ]R2 X 5*^, the leading correction is non-analytic in 
the 't Hooft coupling A, and is approximately of order log(A)^. The free 
energy is expressed in terms of the scalar thermal mass m, which is generated 
by screening effects. We show that this mass vanishes to 1-loop order. We 
then go on to 2-loop order where we find a finite and positive mass squared 
m?. We discuss differences in the calculation between Coulomb and Lorentz 
gauge. Our results indicate that the free energy is a monotonic function in A 
which interpolates smoothly to the N'^/'^ behaviour at strong coupling. 
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1 Introduction 

Realizations of the AdS/CFT duality \T\ exist in several different dimensions. Recently, 
great progress has been made on understanding the gauge theory side of the AdS^/CFT^ 
version of this duality. The original upsurge in interest was generated by the work by Bagger, 
Lambert and Gustavsson now known as the BLG theory. The search for a description 
of the field theory side of the duality, i.e. a world-volume theory on M2-branes has a long 
history, and various obstacles had to be overcome. For example, Schwarz had shown that it 
was impossible to preserve the right symmetries in Chern-Simons theories based on U{N) 
gauge groups [3j . However, Bagger and Lambert managed to write down a Lagrangian with 
all the right symmetries: superconformal symmetry 05*^(814) and parity invariance. In the 
original formulation, it was a Chern-Simons theory based on an algebraic construct known as 
a "3-algebra" , later reformulated as an ordinary quiver gauge theory by van Raamsdonk [1] 
(a possibility which was already discussed by Bandres et al. [5]). However, this theory, 
impressive as it was, turned out not to describe more than (at most) two M2-branes In 
response, various attempts were made to generalize the BLG theory [7]. 

This background set the stage for the breakthrough paper [S] , which successfully general- 
ized the set-up to an arbitrary number of M2-branes. This theory, now known as the ABJM 
model, is an A/" = 6 superconformal Chern-Simons field theory. It builds on previous work 
on superconformal Chern-Simons-matter models [9iLl0j, but with enlarged supersymmetry. 
Crucially, the superpotential in this formulation, while reducing to the BLG superpotential 
for gauge groups of rank 2, did not suffer from the same obstructions as the one used to 
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define the BLG model. In particular, the ABJM model allows arbitrary ranks of the 
gauge group, circumventing the severe gauge group restrictions that unitarity placed on the 
BLG model [11] • The paper [8] includes an analysis of the moduli space for U (N) x U{N) 
gauge groups, showing that M2-branes sitting on a space with a singularity is in fact 
a consistent dual interpretation of the theory. Various checks of this interpretation have 
been carried out. The conjectured JV = 6 supersymmetry of the model was promptly con- 
firmed [HITS]. Progress has been made on understanding the role of monopole operators 
and their relation to the expected supersymmetry enhancement to N' = 8 for Chern-Simons 
levels k = 1,2 [13]. Calculations of the superconformal index match between the strong and 
weak coupling regions [12]. Relations back to the original BLG theory have been established 
in some cases [16]. There is even mounting evidence that the ABJM model is integrable in 
the planar limit [T7HT9] . Chern-Simons theories often arise in condensed matter systems, 
and possible applications include recent studies of the integer and fractional quantum Hall 
effect and Hall transitions [2n] and superconducting M2-branes f2T\ . 

In this note we will study the ABJM model on x at weak coupling and finite 
temperature. Hence, we begin by introducing the ABJM model at finite temperatature in 
section [2j In the dual picture, finite temperature corresponds to a black hole geometry [22] . 
Thermodynamic properties of the theory are most succinctly captured by the free energy. 
This object has been extensively studied in the d = 4, M = 4 super Yang-Mills setting, 
both at weak [23] and strong [231 ES] coupling, including corrections in inverse powers of 
the coupling [25]. At weak coupling, the next to leading order correction is due to thermal 
screening of gauge fields and scalars f27]. One interesting question that we aim to study is 
whether the ABJM theory also exhibits thermal screening. 

In the ABJM case, the strong coupling answer was calculated by Klebanov and Tseytlin 
[25] . and a discussion of corrections can be found in [2B]- On the gauge side, the free field 
theory result was obtained in [8]. Comparing the results at strong and weak coupling, 
there are intriguing differences between the JV = 4 super Yang-Mills theory and the ABJM 
theory. In the ABJM case the gravity result is proportional to N^^'^, which is an artifact 
of M2-branes, while the field theory side has degrees of freedom. Moreover, the entropy 
goes to a constant in the strong coupling limit in the M = 4 super Yang-Mills case, which is 
in stark contrast to the A~^/^ behaviour of the ABJM case. 

One motivation for us is to make progress on understanding how these behaviours come 
about in ABJM theory. Correspondingly, we compute the free energy, including the first non- 
vanishing quantum correction, in section [3l In section [3?T| we show that a naive perturbation 
expansion does not work, which is closely related to the fact that the first non-vanishing 
quantum correction is indeed due to thermal screening of scalars, as we discuss in section 
13.21 Instead, we need to use a method based on resummation of ring diagrams, as explained 
in section 13.31 We then combine our results to write down the full free energy, including 
the leading correction, in section 13.41 Our results are consistent with a free energy which 
interpolates smoothly to the A^^/^ behaviour at strong coupling. We discuss this and other 
issues in section |3l which also contains a short summary of our results. 

This note contains three appendices. In appendix [XJ we list propagators and Feynman 
rules, and explain our notation. Appendix [B] contains the details of the calculation of 
the thermal mass for the scalars. This calculation is a crucial part of this note, but has 
been moved to the appendix due to its lengthy and technical nature. The main gauge 
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used throughout this note is Coulomb gauge, which is often very convenient for thermal 
calculations. Appendix O illustrates some of the difficulties with Lorentz gauge, in particular 
the apparent non-existence of an IR regulating mass for the gluons, which are propagating 
in Lorentz gauge. 



2 ABJM at finite temperature 

The ABJM model is a three-dimensional Chern-Simons-matter theory with gauge groups 
U{N) X U{N) and SU{N) x SU{N), with J\f = 6 superconformal symmetrjl^. Other choices 
of gauge groups exist, most notably those which are dual to orientifolded geometries [30]. 
At large N, the U{N) x U{N) theory is believed to be dual to M-theory on AdS^ x S'^/Zk. 
Henceforth, we will restrict attention to the U{N) x U{N) theory. In this section, we will 
also define the ABJM model at finite temperature. 

The ABJM model is defined by the following actior|§. 



S 



k 

2^ 



d X 



e'^HT{-^AidjAk + \AiAjAk + '^-A.djAk - \A,AjAk) 



(2.1) 



We are using the notation of [12], but we have rescaled all matter fields by Y"^ — t- y -^Y"^ and 

ipA to be able to factor out the Chern-Simons level k as an overall normalization 

of the action. In addition, we have performed a Wick rotation to Euclidean space. The 
Dirac matrices are (7*)^'' = (— cr^, cr^, o"^), where a* are the Pauli spin matrices. 

The action f l2.ip exhibits an A/" = 6 supersymmetry enhancement, owing to an underlying 
SU{4) R-symmetry which acts on the matter fields Y"^ and ipA = 1,2,3,4). As shown 
in |8l|T2], this symmetry can be made manifest in the bosonic and fermionic potentials by 
writing them on the form 



V 



bos 



tr 



yIy^yIy'^y^y^ 



+ ay'^yIy^yIy^y}; 



QY'^Y]jY"YlY'-Yc 



V 



fcrm 



i tr 



'\yU^^ 

- e^^^^Yli^sYUn + eABCDY^^^^Y^^^'' 



Y),Y''^'^^B - Y^YlijB^^'' + 2F^FiV^^V^ti? - 2rjF^^tA^^ 



(2.2) 



(2.3) 



The matter transforms in the bifundamental representation of the gauge group. Hence, the 
covariant derivative is ViY^ = diY^ + iAiY^ — iY^Ai. The gauge fields (henceforth referred 
to as gluons) Ai and Aj {i = 1, 2, 3) are in the adjoint of the left and right U{N), respectively. 
The field content of the ABJM model is summarized in figure [H 



^An introduction to the ABJM model is given in 

^ For a discussion on pure spinor superfield formulations, see |31] . 
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The action (12. ip has a U{N) x U{N) gauge-invariance, and needs to be gauge-fixed. We 
gauge-fix by adding 



J £x tr (^('9;.^'^)' + - (c d.V' c + d d.V' c) 



(2.4) 



to the action, S ^ S + Sg. Two possible gauge choices are Lorentz gauge (/i = 1, 2, 3) and 
Coulomb gauge (/i = 1, 2). We typically use Coulomb gauge in the main text, and comment 
on the use of Lorentz gauge in footnotes. ^ and ^ are dimension one gauge-fixing parameters 
for A and A, and c and c are the corresponding ghosts. Hence, the covariant derivatives are 
^^^lC = d^c + iA^c - icA^, and P^c = d^c + iA^c - icA^. 

To study the thermal theory, we compactify the time direction on a circle of length 
L = 1/T, where T is the temperature. The action (12. ip and (12.40 then defines the theory on 

X S*^. Finite temperature breaks both supersymmetry and conformal invariance. A generic 
field has boundary condition x^, + L) = (— x^, x^) along the compactified 
direction. In particular, the scalars, gluons and ghosts have periodic boundary conditions 
(z/ = 0), and the fermions are antiperiodic {ly = 1/2), i.e. 

X^ J X^ J X^ ) • 

Ai{x\x^,x^ + L) = +Ai{x\x'^,x^), 
c(x\ x^, x^ + L) = +c(x\ x^, x^). 

Propagators and vertices can be derived from the action (12. ip and (12.41) . These Feynman 
rules are collected in appendix \M This appendix also contains an explanation of some of 
our notation. 



3 Perturbative free energy 

We now want to compute the free energy, including the leading correction, and analyze its 
coupling dependence. We will show that to obtain a finite answer, we will need to reorganize 
perturbation theory, taking thermal screening of the scalars into account. 
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3.1 Naive perturbation expansion 

In this section, we will compute the free energy in a naive perturbation expansion. We will 
show that perturbation theory breaks down between two and three loops. 

Assuming that the free energy is analytic in the 't Hooft coupling A = we write 

F{\) = Fi + XF2 + X^F3 + --- . (3.1) 

We begin at one-loop order, using Coulomb gauge. The free energy receives contributions 
from scalars and fermions (see figure [2]) , which sum up to 

F, = (sUo - sUy^) = -iV^T^^. (3.2) 



where 



n=£Z+v 

A, = -T^^, (3-3) 
, T.33C(3) 

corresponding to contributions from fields with periodic (z/ = 0) or antiperiodic {v = 1/2) 
boundary conditions. Notice in (13. 2p that the matter fields each have eight real degrees of 
freedom. Relative signs arise for the anticommuting fermions. To verify that gluons and 
ghosts do not contribute, note that 

det^/2(±6'^-^9. + ^94) = -|=. (3.4) 

An arrow denotes that the object only has components along the 1,2 directions, i.e. no 
dependence on p^. Thus, this contribution vanishes in zeta function regularization. It is 
the Chern-Simons terms which make it possible to find a gauge in which gluons and ghosts 
are non-propagating. In the zero-temperature limit T — t- 0, the free energy goes to zero 
(independently of gauge choice). At zero temperature (i.e. without compactifying the time 
direction), contributions from bosonic and fermionic degrees of freedom cancel automatically, 
as is required by supersymmetry. The result (13.21) was already derived in [S]. 



^ In Lorentz gauge, there will also be contributions from gluons and ghosts, since 

92 



dei'/\±e^^^du + ^a.a,) = -j=, (3.5) 



becomes 



9/11 \ 9 •^7C(3) 



Fi=N'[ 8-Ao - 8-A1/2 + 2Ao - 2Ao ) = ^N'T^^^. (3.6) 



Notice that the two types of gluons and ghosts (which are anticommuting, hence the relative sign) now do 
not vanish individually (as in Coulomb gauge), but instead cancel each other. In neither gauge is there a 
net contribution to the free energy from gluons and ghosts. 
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Figure 2: One-loop contributions to the free energy arise from (a) scalars and (b) fermions (and in 
Lorentz gauge, also from (c) gluons and (d) ghosts). 




+ 






Figure 3: Two-loop contributions to the free energy. 



Using the interactions f lA.ip and propagators (1A.4I) . we can now try to find the perturba- 



tive corrections to the free energy. At two-loop order, the free energy receives contributions 
from the connected and one-particle- irreducible diagrams shown in figure [31 All the diagrams 
vanish individually, due to combinatorics. This is true both at zero temperature and finite 
temperature, and for both Coulomb and Lorentz gauge. Hence, 

F2 = 0. (3.7) 

In view of fl3.7p . we might think that the first non- vanishing correction terms to the free 
energy are of order 0{X^). However, this conclusion assumes that the perturbation expansion 
(13. ip is well-defined around A = 0, and this is not true. The obstruction consists of infrared 
divergences, appearing at three-loop order. For example, the diagram shown in figure H] is 
proportional to 

The infrared divergences do not cancel. One possible way to regularize the theory is to go 
to finite volume. However, we are expressly interested in the non-compact theory, where the 
divergences should instead be cured by summing over ring diagrams. We will explain how 
this works in the following sections. 



3.2 Thermal screening 

In section 13. we found that perturbation theory breaks down between two and three loops. 
A similar phenomenon has been analyzed previously in QCD [32HS1] and in super Yang- 
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Figure 4: A three-loop contribution to the free energy. 



Mills [27], where scalars and gluons are screened by quantum effects. The ABJM theory 
exhibits a similar behaviour. 

To proceed, we must reorganize perturbation theory in a way which regularizes the in- 
frared divergences in the propagators [32] • In momentum space, the quadratic part of the 
action (O) and (EID is 



^2 = — 

1 

2? 



(2^ 



+ + + c{p)p^p''c{-p) + c{p)p^p^'c{-p) (39) 

+ Ylip)p'Y^i-p)]-^ f^T J2 tr[p^^ip)fM-p) 



We now write the action as S* = (5* + 55*2) — 55*2, where 



5S, = ^ 



(Pp 



277 J (27r)^ 



tr 



^Wl{p)ml.Y\-p) 



(3.10) 



We henceforth treat —882 as a perturbation to {S + 5S2)- Effectively, this means that 
infrared divergences are regularized by the generation of a thermal mass my in the propagator 
for the scalaro 

1 

(3.11) 



2tt 

(Ynip) YU-P)) = — 



kT fp + C0I + my 

The generation of thermal masses is similar to the renormalization of the electric charge in 
QED, where virtual electron-positron pairs are created and cause a dielectric screening effect 
from the vacuum. A thermal mass m corresponds to generation of a finite static screening 
length r = ^ [33]. 

Technically, the scalar thermal mass is obtained by computing the one-particle- irreducible 
contributions to the self-energy in the static limit (the limit of no external momentum). The 



In other gauges, thermal masses may be generated for all fields with bosonic boundary conditions. 
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details of the calculation (in Coulomb gauge) are collected in appendix |Bl The final result 
is given by equations (]B.13|) and (IB.lSp . quoted here for convenience: 

ml{\) = (27rr)V(A), 
A^) = log(/.)^ + 0{X' log(A)). ^^-''^ 

A few comments are in order. Since the first non-vanishing corrections (i.e. at two loops) 
to the thermal mass are themselves infrared divergent, the thermal mass must be computed 
self-consistently by requiring that successive higher order perturbative corrections do not 
shift the pole in the propagator (cf. D'Hoker's original calculation for QCD3 [35j). Thus, the 
diagrams will themselves depend on the regulator that we are trying to compute. Notice that 
the right-hand side of equation fl3.12p depends on /i^(A). However, to a good approximation 
(for very small A), we can replace log(/i) by log(A). Notice that my is then approximately 
of order A^ log(A)^, rather than just A^. 

In principle, it would have been useful to verify that we are indeed expanding around 
the correct vacuum of the finite temperature theorjU. Here, we will be content with the 
observation that the mass squared my is manifestly positive. It would also have been useful 
to verify our results by an analogous treatment in Lorentz gauge. In Lorentz gauge, gluons 
are propagating, and would apparently need a regulator to yield a finite answer. However, no 
such IR regulating mass exists for the gluons up to two-loop order, as we show in appendix 

o 



3.3 Reorganized perturbation expansion 

We have concluded that finding the leading correction to the free energy requires reorganizing 
perturbation theory. Thus, we need to replace (13. ip by 

F(A) = Fi(A) + AF2(A) + A^F3(A) ■ • • , (3.13) 

where -Fi(A), F2{\), F-^^X), ... are now computed using the renormahzed scalar propagator 
(13. lip . There is a A dependence in -Fi(A), -F2(A), F^i^X), .... The reason is that there is a A 
dependence in the thermal mass f l3.12p . In fact, it is this A dependence which gives rise to 
the first non- vanishing (and non-analytic) correction to the free energy, as we will now show. 

We begin with the resummed one-loop contribution -Fi(A) (see figure [5]). Still working in 
Coulomb gauge, we find 

i^i(A) = (^S^Aoim'y) - 8^Ai/2(m|)^ , (3.14) 

which generalizes the free field theory result (13.20 to the case of non- zero masses. Fields with 
periodic boundary conditions {u = 0) and mass m contribute 

Ao{m')=Tj2 I 7^1og(p-^ + (27rrn)2 + m^)= (3.15) 

= [-4C(3) - m'T-^ log (m'T-2) + m^T"' + O (m^T"^)] , 
^ I thank E. Witten for drawing my attention to this issue. 
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Figure 5: Free energy ring diagrams, 
and fields with antiperiodic boundary conditions {u = 1/2) contribute 

Ay.im') = T J2 + {2nTnf + m') = (3.16) 

= ^ [3C(3) - m^T-^ log(4) + O (m^T"^)] . 

In the free field theory limit A — t- 0, we recover the leading order results (13. 2p and (13.31) . i.e. 
linim-j^o ^o("^) = ^0 and linim^o ^i/2("^) = ^i/2- Note also that there is no zero mode for 
fields with antiperiodic boundary conditions, so no thermal mass will be generated for them. 
The dependence on in (I3.16P is analytic. Moreover, it follows from (I3.15P and (I3.16P that 
the leading order correction will anyway come from fields with periodic boundary conditions 
(i.e. the scalars). 

The leading correction to the free energy is fully contained in Fi{X). Neither -F2(A) nor 
F^IX) contributes, for the following reasons. F2{X) still receives contributions from the same 
diagrams as before, shown in figure [31 Hence, it still vanishes for combinatorical reasons. On 
the other hand, -F3(A) is not expected to vanish. However, all contributions to the leading 
correction are already included in -Fi(A), due to the definition of the thermal mass. Notice 
that we obtain 3-loop diagrams by closing the scalar propagators in the diagrams in figures 
[HI [^ and [TUl Any such prospective additional contributions to -F3(A) will be cancelled by 
counterterms. Equivalently, D'Hoker's self-consistent treatment requires that the pole in the 
scalar propagator is not shifted by higher order corrections [35] . 

3.4 Dependence on the coupling 

We have argued that the leading correction to the free energy is non-analytic in the coupling, 
due to thermal screening of scalars. We now combine these results to analyze the coupling 
dependence of the free energy. 

The free energy density, including the leading correction, is 

F = -N^T^fiX), (3.17) 
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Figure 6: Free energy and scalar thermal mass. 



7rT2 



2^2 



(3.18) 



as follows from combining equations f l3.13p . f l3.14p . f l3.15p and f l3.16p . and with the scalar 
thermal mass my(A) defined by equation f l3.12p . Using the approximation log(/i) ~ log(A), 
we get the convenient closed form answer 



/approx(A) 



^ + f^A^log(A)^ + 0(A^log(A)^ 



(3.19) 



We observe that the leading order correction to the free energy is (approximately) of order 
A^log(A)^. Such a behaviour has also been observed in other three-dimensional Chern- 
Simons- matter theories with lower supersymmetry [9]. 

Not using the approximation log(yu) ~ log(A), equation f l3.12p can be solved numerically. 
The numerical solution, and the resulting free energy, are shown in table [T] and plotted in 
figure El The solutions are not plotted beyond A = exp(— 7) ~ 0.001, since the next to 
leading order correction will be negligible only so long as log(A) is sufficiently large. The 
free field theory result for the free energy is /(O) ~ 2.6784. 

We see that our results indicate that the free energy is a monotonic function in A which 
interpolates smoothly to the N'^/'^ behaviour at strong coupling, where the free energy is 
given b}@ ^ 



/s(A) 



2^/2 2 1 

9 



(3.20) 



For comparison, in A/" = 4 SYM the analogous results in the weak and strong coupling regions are 



Fsym(A) = -—N^V^T^giX), 



5wcak('^) 
^strong (A) 



3 45 1 

4 + 32^^'^^a57^ 
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9 fi7797 


nil c^8/i 


0.0007 


2.67698 


0.15049 


0.0008 


2.67668 


0.18864 


0.0009 


2.67636 


0.23011 


0.0010 


2.67603 


0.27473 



Table 1: Free energy and scalar thermal mass. 

4 Discussion 

In this note, we have computed the free energy in ABJM theory on M^x S\ given in equations 
fl3.17p and fl3.18p . The free energy is expressed in terms of the scalar thermal mass fl3.12p . 
whose generation is due to screening effects. The numerical solution is shown in table [1] and 
plotted in figure O Our answer for the free energy includes the first non- vanishing quantum 
correction, which is approximately of order A^log(A)^ (see equation fl3.19p ). The reason for 
the non-analytical dependence on the coupling is that the IR divergences had to be cured 
by using a technique based on resummation of ring diagrams. 

Interesting differences appear compared to the calculation in A/" = 4 super- Yang Mills 
theory (SYM) j25l[27] . One important difference is that in SYM, the first non- vanishing 
correction is analytic in the coupling, and only after that do the non-analytic contributions 
appear. Another difference is that the self-consistent treatment is not necessary, since the 
1-loop diagrams which define the thermal mass are regular. For us, already the first non- 
vanishing correction is non-analytic in the coupling, and the self-consistent treatment is 
required (note that the right-hand side of (13.121) explicitly depends on the mass itself). 

It would be interesting to compute higher order corrections to the free energy. However, 
even in the reorganized theory, the perturbation expansion is eventually expected to break 
down (at finite temperature) [36]. A related observation is that there are some indications 
[311137] that even the self-consistent treatment proposed by D'Hoker [35], which we follow, 
does not always eliminate all infrared divergences, due to differences between electric and 
magnetic masses. Our calculation of the free energy shows explicitly that the reorganized 
theory furnishes a finite answer to order A^log(A)^, but there is no guarantee that this 
finiteness persists to higher orders. 

Another interesting question is the issue of phase transitions. In SYM, there were some 
initial signs that the system may go through a phase transition between strong and weak 
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coupling [38]. However, the present understanding is that the free energy is a smooth and 
monotonic function for all values of A [39], and more recent evidence also exists [IQ]. However, 
even if the phase diagram on flat space is trivial in SYM, there are examples of three- 
dimensional field theories with possibly different behaviour, e.g. QCDg [35] • Nevertheless, 
we would expect ABJM to behave similarly to SYM in this respect. Our results for the free 
energy and the scalar thermal mass are only valid so long as log(A) is sufficiently large, but 
they still allow us to conjecture that /(A) is indeed smooth and monotonically decreasing 
between the weak and strong coupling regions also in ABJM theory. 

A non-trivial phase space behaviour is expected on compact spaces. In the gravitational 
black hole description, the existence of a deconfinement (Hagedorn) transition has been 
established [22]. SYM theory on a compact space was analyzed at zero 't Hooft coupling by 
Sundborg by counting gauge-invariant states on a 3-sphere using Polya theory ^1]. These 
results were verified and extended by Aharony et al. [12], and later complemented by one- 
loop corrections [13]. Similarly, using Polya counting, Nishioka and Takayanagi found a 
deconfinement transition in ABJM theory compactified on a 2-sphere [T^, and further study 
of such questions in ABJM would also be interesting. 
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A Propagators and vertices 

In this section, we list propagators and interaction vertices, which can be derived from the 
action and M . 
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Figure 7: Interaction vertices. 



The interaction vertices are as follow^ ( see figure [7]) . 



1 

— tr 
3 



V 



tr r^iM^rj + yIa'A'y^ - ia'yIa'y^ , 

ztr [y-jy^^t^^/'B - r^FjVsV'^^ + TY^yI^^^^'^^ - 2Y)^^^\}'^^^\}b 



T/gh 



Igiifc^r 
3 



% tr 
tr 
z tr 



AiAjAj^ — AiAjA)^ 
d,Y\{A'Y^ - Y^A') + diY^{A'Yl - 

a^c(A^c - cA'^) + 9^a(>c - cA^) . (Coulomb gauge: = 1, 2) 



(A.l) 



Using the Fourier transform 



d X-, 



c(s 



\x\ 



\p\ 



3-s 



(A.2) 



c s 



"KS 

T 



47rr(2 -s) sin f- 

we can derive the momentum space propagators (or read them off directly from (13.91) ) in 



^ In Lorentz gauge, ^=1,2,3. 
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Coulomb gaugqj: 





-P) 




-P) 




-V) 




-V) 


(Cnip) C_„( 


-V) 


{Cn{p) C-n{ 


-P) 



2tx 


1 




kf 






2n 


f ^ 
■I a 




kf 






27r 
kf 




1 


2n 
kf 




1 


2n 


1 




kf 
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1 




kf 


p2 





(A.4) 



We use the Landau choice C, = ^ = 0. Deha functions in colour and flavour indices are 
suppressed. In writing down the propagators, we use the notation that p = {pi,P2), {P3)n = 
Un = 2'KTn. The modulus is ^/p^, where p^ = p^ + (a;„)^. The index n refers to different 
Fourier modes along the compactifled direction. The mode expansion must adhere to the 
boundary conditions (12. 5p . Hence, n G Z + z/, where v is 1/2 for the fermions ipAi and 
otherwise zero. For example, the scalars and the fermions are expanded as 



We also make use of the symbols 



(A.5) 



A{p) 

Mp) 



1 



p'^ + UJ^ + my 



2 ' 



(A.6) 



Aa{p) = AM = ^, 



^ In Lorentz gauge, the propagators for gluons and ghosts are 

2n 



kl p + ui„ 



{Cn{p) C-n{-p)) = - 



kT jp + 
2n 1 



(A.3) 



kT + a;2 



{Cn{p) C-n{-p)) = ■ 



kT jp + uj^ 
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(where p = {p,27iTn)), corresponding to the generic, scalar, gluon and ghost propagator, 
respectively. Notice that As{p) includes the mass renormalization, as in (]B.14p . The symbols 



and /\c(p) are equal, but it is still useful to remember their respective origins. We use 
the subscript / to denote fermionic boundary conditions. Explicitly, Pf = {p, 27rT(n + 1/2)). 
For example. 



(pi+P2)^+W^^+n2+l 

An arrow on top of an argument denotes setting p^ — > 0. For example. 



Aip) = -. (A.8) 



B Scalar thermal mass 

In this section, we compute the scalar thermal mass to leading order in Coulomb gauge. 
Stated more precisely, our objective is to compute the scalar self-energy in the limit of 
vanishing external momentum (the static limit), defined by the sum of all one-particle- 
irreducible diagrams. 



-O = -T lim {p,0) = -27r/cT ( A) . (B. 1) 



We will find that the one-loop answer vanishes. We therefore need to proceed to two loops 
to find a finite answer. 

We begin with the one-loop calculation. We will organize the computation in the following 
way. There are two distinct types of diagrams which contribute at one-loop order. These 
diagrams contain interaction terms of the types 

type (4): {Vi + V^) , 

type (33): [V^ + + Vf + Vf)', 



in the notation of flA.l|) . There are 9 distinct terms of type (4) and 78 distinct terms of type 



(33). From flA.ip and (IA.4|) . it follows that in the static limit (no external momentum), all 



one-particle irreducible diagrams vanish: 
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?^(70/ = 0, (B.2) 



~ e'^^5^j = 0, (B.3) 

(2) ~ e^^'^g.g, = 0. (B.4) 

Hence, we need to proceed to two-loop order. At this order, four different types of 
diagrams contribute, due to interaction terms 

type (6): Ve, (B.5) 

type (44): {Vi + Vl)\ (B.6) 

type (433): (V? + Vl) (v,^ + + Vf + Vf)', (B.7) 

type (3333): {v^ + + Vf + V^y, (B.8) 



still written in the notation of flA.ip . There are 4, 45, 702 and 1365 distinct terms of type 
(6), (44), (433) and (3333), respectively. 

We first consider a subset of the (433) diagrams. Thermal ABJM theory is severely 
plagued by IR divergences. However, there are some cancellations. In particular. 
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4 

Igiqu 92) = - 2 (gi ■ g2)^^A(gi)^A(g2)^A(gi + 52)^ 

+ ^g?g^Z\^(gi)Z\^(g-'2)Z\^(g-'i + q^Y 
+ ^{qi ■ q2fA,{qi)A,{q2)AA{qi + 

- ^gt^t^clgO/^clgD/^Algi + gl)' = o, 



in the notation of flA.6p and (lA.Sp . Hence, all ghost contributions cancel already in the 



integrand by equal but opposite gluon contributions. Another way to see this is to notice that 
(IB.9P contains a piece which is constructed from 1-loop corrections to the gluon propagator. 
This piece cancels by itseli@: 




We now consider all remaining diagrams. One-particle-irreducible and non-vanishing dia- 
grams contributing in the static limit are shown in figures [HI ID and [TDl Adding together all 
contributions, we fincrl 



where 



- l™/''«0) = 2"^= / Y: j Y: Hiu,.), (B.IO) 

Hqu 92) = hiqu 12) + h4iqi, q2) + /433(gi, ^2) + /ssssl^i, ^2), 

/6(gi,g2) = -4:2As{qi)As{q2), 

/44(gi, q2) = -32(iAsiqi)AAiq2f ^^'^^^ 
+ A{qi ■ q2)As{qi)AA{q2)AA{qi + 52) 
+ A^^As{qi)AA{q2)AA{qi + ^2) 
- 112(g/i ■ qf2)A{qfi)A{qj2)As{qfi + g/2), 



The same result holds in Lorentz gauge. 

In the flat case and for general momenta, the scalar self-energy to two loops was calculated in [T5] . 
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Figure 8: These type (6) and (44) diagrams are non-vanishing in both Coulomb and Lorentz gauge. 




Figure 9: These type (433) diagrams are non-vanishing in both Coulomb and Lorentz gauge. 



18 



32 

^433(gi, ^2) = - y (gl ■ g2)^^(g/i)^(g/2)^A(gi + 

64 

- y'?i('?i ■ ^2)^(g/i)^(g/2)^A(gi + ^2)^ 

32 

- —(f,qlA{qfi)A{qf2)Mqi + ^2? 
2 

- 3(91 ■ ?2)(gi ■ q2)AA{qi)AA{q2)As{qi)As{q2) 
2 

+ -(gi ■ q2fAA{qi)AAiq2)Asiqi)Asiq2) 
16 

+ yg2?1Z\^(g2)M5(gi)Z^5(g2) 
16 

- —(f2f^AA{q2fAs{qi)As{q2) 
4 

- 2(91 ■ 92)(gi ■ g2)'^A(g2)^A(gi + q2)As{qi)As{q2) 
4 

- ^(^L ■ q2YAAiq2)AA{qi + q2)As{qi)Asiq2) 
4 

- 2^2(^1 ■ q2)AAiq2)AAiqi + g2)^5(gi)^5(g2) 
4 

- ■ g2)^A(g2)^A(gi + q2)As{qi)As{q2) 
4 

- 3?2^^A(g2)^A(gi + q2)Asiqi)Asiq2) 
4 

+ 2g^'?^^A(g2)^A(gi + q2)As{qi)As{q2) 
32 

+ y (gi ■ ?2)(gi • g2)/\A(gi + q2fAs{qi)As{q2) 
32 

+ y^i(^i ■ ^2)^A{qi + q2fAs{qi)Asiq2) 
32 

+ y'?i(^i ■ ?2)^A(gi + g2)M5(gi)Z\5(g2) 
16 

+ yg2g-^Z\^(gi + g2)M5(gi)Z\5(g2) 

16 

- —q^^q^^AAiqi + g'2)M5(gi)Z\5fe) 
16 

+ yg2g-^Z\^(g, + g2)M5(gi)Z\5(g2) 

16 

+ yg?g1Z\^(gi + g2)M5(gi)Z\5(g2), 

in the notation of (1A.6I) . (1A.7I) and (1A.8I) . These expressions have a tendency to become 
more and more involved as the number of vertices increases. We omit spelhng out the full 
expression for /3333(gi, q2)- As we will see below, it has no zero modes and will not be of any 
further interest to us. 
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Figure 10: These type (433) diagrams vanish in Lorentz gauge (in the static hmit) but not in 
Coulomb gauge. 



We expect the leading order result to arise from the most IR divergent piece. This piece 
(specifically, the piece proportional to A^log(/x)^) can be extracted from the zero modes 
(among which many terms cancel). The zero modes will provide reliable results to leading 
order in A, but will also contribute to subleading terms. However, sub leading terms can 
and will receive contributions from other modes. Hence, we will only collect the leading 
order result. We also take into account that g! and q2 are dummy variables that can be 
interchanged in individual terms. Terms which do not depend on both ni and 77-2 vanish in 
zeta function regularization. Taking all this into account, the integrands become 

Iei(lu<l2) = -A2As{qi)As{q2), 

(12) = ^ - 2g?g1 - (gi + q2)\qi + gl)'] AA{qi + g2)M(g->i)Z\(g->2) + 



^3333(^1.^2) 



+ 7T [4gt(g'i + q2f - 2g?gt - (^i + <l2f{qi + gl)'] ^A{qi + q2f As{qi)As{q2] 



3 
0, 



still written in the notation of (IA.6p . flA.7p and flA.Sp . We still need to do the integrals 
over qi and g2- The theory can be regularized in the ultraviolet (finite temperature only 
changes the theory at scales > j^). Thus, the leading order result can be collected from the 
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IR divergent end of the integrals. The result is 

- hmM'ip, 0) = 2nkT'\' J J [/"(gl, q,) + I^^^^(q-[, q,)] + higher modes 



27rkT^X^ 



-42 ^og{,f + f ^ loM^ - log(^)^ 



0(AMog(A)) 



118 



log(/i)^ 



_3(27r)2 

Comparing to fIB.ip . we identify 



+ 0{\Hog{\)). 



3(277)^ 



(B.12) 



(B.13) 



Notice that the right-hand side itself depends on /U^(A). This is due to the fact that consis- 
tency (IR finiteness) requires that we use renormalized propagators (1B.14P when we compute 
the values of the self-energy diagrams Summing up the one-particle-irreducible contri- 
butions in the conventional geometric series 



+ •- 



we find the renormalized scalar two-point function (external lines are included) 



27r 



kT f + ujl + ml 



2 ' 



(B.14) 



where the scalar thermal mass is 



m^(A) = (27rT)V(A) 



and /i^(A) is given by equation flB.13p . 



(B.15) 



C Gluons 

In this section, we compute the static self-energy of gluons and ghosts in Lorentz gaugi^^^' 
We will find that to two-loop order, no thermal mass is generated. 



Note that Coulomb gauge is our main gauge choice throughout ah other parts of this note. This section 
is primarily meant to illustrate the problems with Lorentz gauge, in particular the non-existence of an IR 
regulating mass in a gauge which leads to propagating gluons and ghosts. 
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We first consider one- loop diagrams. Using (lA.ip and (lA.4p . we find that the following 
one-particle- irreducible diagrams contribute to the gluon self-energy in the static limit. 




(C.l) 




^^kTX I ^jY.^,q,As{q)\ (C.2) 



(C.3) 




+ 



(C.4) 



written using the definitions in (lA.6p . The sum of gluon and ghost loops cancel exactly if 
we use uncorrected propagators. In principle, the result flC.4p may receive corrections from 
renormalized propagators [35]. It turns out that diagrams contributing to renormalization 
of A and A are equal. We can sum up fIC.ip . flC.2p . (IC.SP and flC.4p either by doing the sum 
or the integral first. If we do the sum first, the summation formulas 



and 



y— 



TT coth(7rZ\) 
1 TT tanh(7rz^\) 



n2 + Z\2 



A 



E 
E 



'n2 + Z\2 



7rcoth(7rZ\) 7rcoth(7rZ\) 

"2 ~ ^P,0 TTT^ 1" 



2Z\3 



2Z\ 



, (n2 + A^y 



7rtanh(7rZ\) 7rtanh(7rz^) 



2Z\3 



2Z\ 



(C.5) 



(C.6) 



(for p G {0, 1,2}) are useful. The formulas fIC.Sp and flC.6p can be derived using a method 
similar to the one in [H], by considering the contour integrals 



= It + Jb + /c 

= It + 4 + /c 



dz 

—f{z)'KCot{nz), 
Zm 

dz 
Ziri 



(C.7) 
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c 



At 



Figure 11: Integration contours and poles. 



Poles and integration contours t, b and c are shown in figure [TT] (for tlie case witli a pole at 
the origin). Note that the function it cot{7rz) has simple poles at z = n, n & Z, and that the 
function — 7rtan(7rz) has simple poles at z = n, n E Z + ^. The summation formulas fIC.Sp 
and flC.6p then follow by setting f{z) gi{z) and f{z) — (72(2), respectively, in flC.7p . The 



functions gi{z) and g2{z) are defined by 

where p G {0, 1, 2}. The functions gi{z) and g2{z) have simple and double poles, respectively, 
at z = ±iA. 

Either way, dUl]), (ES), (El and ([03]) sum up tc@ 




= limi7,,(p,0) = 5i,35j,3kT^n{X,fj.), (C.IO) 



where i7(A,/i) = |^ log(/i) (to lowest order in A). Hence, the gluon two-point function 
I thank J. Maldacena for discussions on the expected form of the gluon self-energy. 
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becomes (upper sign corresponds to Ai, lower sign to Ai) 



kT 

— ( ±eijkPk + i^PiPj ) + Si36j3kT'^n 



In 



(2 



71 



{kTf 

where the tree-level propagator is 



2i 



kT'^n 



{p'y 



-P1P2 Pi 

\ 



P1P2 

2 




0/ 



kT\-^ 
2^) 



( 



±UjkPk + ^PiPj 



2tt 1 
kT p'^ 





P3 



~P-i P2 

-Pi 



(C.ll) 



(C.12) 



\-p2 Pi / 

and the arrow denotes the operation of inverting the operator and then using the Landau 
gauge choic^ = 0. 

Comparing fIC.lip to flC.12p . we note that no mass is generated and the correction is 
not sufficient to suppress the IR divergences. Hence, we need to proceed to higher orders 
to find an IR regulator for the gluon, if it exists. At two-loop order, only two diagrams can 
contribute for combinatorical reasons. 




= 27rkT\^f [{q2)iejki + (g2),Qfcz] iqi)kiq2)i^A{qi)MQ2?MQi + Q2), (C.13) 

where qi = (gi,27rTni) and q2 = (g2,27rTn2). The last equality sign comes from doing the 
loop integrals. Hence, the static self-energy vanishes at two-loop order. 

Ghost self-energy diagrams always contain the vertex structure Vf^ (as defined in (lA.ip ). 
Using this fact and the static limit, it can be verified that to two-loop order, no thermal 
mass is generated for the ghosts, 

ml = 0. (C.14) 
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